REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden  to  Deparbiient  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis 
Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a 
collection  of  information  if  it  does  not  display  a  currently  valid  0MB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE  3.  DATES  COVERED  (From  -  To) 

10-02-1998  Paper 


4.  TITLE  AND  SUBTITLE  5a.  CONTRACT  NUMBER 


Modeling  of  Interfacial  Fracture  in  Incompressible  Materials 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 
ANDADDRESS(ES) 


5d.  PROJECT  NUMBER 

2302 


5e.  TASK  NUMBER 
M1G2 


5f.  WORK  UNIT  NUMBER 


8.  PERFORMING  ORGANIZATION 
REPORT 


Air  Force  Research  Laboratoiy  (AFMC) 
SPARTA 
5  Pollux  Drive 

Edwards  AFB  CA  93524-7048 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


10.  SPONSOR/MONITOR’S 
ACRONYM(S) 


Air  Force  Research  Laboratory  (AFMC) 
AFRL/PRS 
5  Pollux  Drive 

Edwards  AFB  CA  93524-7048 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 


11.  SPONSOR/MONITOR’S 
NUMBER(S) 

AFRL-PR-ED-TP-1998-04 1 


14.  ABSTRACT 


20020115  079 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

Unclassified 


b.  ABSTRACT 

Unclassified 


c.  THIS  PAGE 

Unclassified 


17.  LIMITATION 

18.  NUMBER 

OF  ABSTRACT 

A 

OF  PAGES 

19a.  NAME  OF  RESPONSIBLE 

PERSON 

C.T.  Liu 


19b.  TELEPHONE  NUMBER 

(include  area  code) 

(661)275-5642 


Standard  Form  298  (Rev.  8-S8) 

Prescribed  by  ANSI  Std.  239.18 


MODELING  OF  INTERFACIAL  FRACTURE  IN  INCOMPRESSIBLE  MATERIALS  WITH  VARYING 

MODULUS  MISMATCH’ 

T.  C.  Miller 
Sparta  Incorporated 

2  Draco  Drive,  Edwards  Air  Force  Base,  CA  93524 

ABSTRACT 

Numerical  modeling  is  used  to  evaluate  the  effects  of  modulus  mismatch  on  interfacial  fracture.  Different 
modulus  ratios  are  considered,  as  are  different  mode  mixities.  The  magnitudes  of  the  complex  stress  intensity 
factors  are  evaluated  using  the  energy  domain  integral  approach,  and  the  phase  angles  are  measured  using 
“olation  of  bond  Infraction  data  to  r  =  0.  The  results  indicate  that  moderate  changes  m  ^e  moddus  ratio 
have  only  a  small  effect  on  either  the  magnitude  or  phase  angle  of  the  complex  stress  intensity  factor 
predictions  confirm  earlier  experimental  conclusions  and  suggest  that  the  elastic  mismatch  is  not  strong 
in  determining  the  interfacial  fracture  of  plane  strain  incompressible  materials. 

INTRODUCTION 

The  use  of  fracture  mechanics  to  assess  and  predict  crack  development  in  solid  propellants  has  led  to  cost 
savings  by  enabling  better  service  life  predictions  for  solid  rocket  motors.  Both  Imear  e  astic 
fracture  mechanics  concepts  have  been  used  to  predict  crack  propagation  m  solid  propellants  A  frequent  site  for 
the  initiation  and  growth  of  cracks  in  solid  rocket  motors  is  the  interface  between  the  propellant  and  the  rubbj 
liner.  TTiis  study  Examines  this  situation,  namely,  that  of  an  interfacial  crack  lymg  between  two  incompressible 
materials  and  subjected  to  plane  strain. 

A  previous  experimental  work  examined  a  similar  situation.'  Both  homogeneous  and  bimaterial  specimens 
were  considered.  In  bimaterial  specimens,  the  elastic  moduli  of  the  two  materials 

Additional  results  also  indicated  that  larger  differences  between  the  two  elastic  moduli  had  little  effect  on 
complex  stress  intensity  factor. 

The  modulus  ratio,  E,/E„  depends  on  the  materials  used  in  the  motor  arid  different  values  need  to  be 
considered.  Typically,  the  stiffness  of  the  propellant  will  exceed  that  of  the  liner  by  a  factor  of  about  3^  If 
can  be  established  that  this  ratio  only  weakly  affects  the  complex  stress  intensity  toor,  then  computations  for 
interfacial  cracking  of  different  motors  can  be  combined  in  a  simplified  analysis.  The  present  work  uses 
numerical  methods  to  vary  E2/E2  from  1  to  6  to  study  this  issue. 

DISCUSSION 

FTET.D  EXPRESSIONS 

Figure  1  shows  a  crack  and  corresponding  coordinate  system  lying  along  an  interface  between  two  distinct 
linear  elastic  isotropic  materials.  If  the  materials  are  both  incompressible  and  plane  strain  conditions  exist,  then 
the  bimaterial  parameter,  s,  normally  associated  with  interfacial  cracks,  vanishes  This  is  a  ^ 

the  more  general  interfacial  cracking  problem  and  it  has  field  expressions  that  closely  resemble  those  of  cracks  in 
homogenous  bodies  subjected  to  mixed  mode  loading.  One  consequence  of  the  vamshmg  of  e  is  that  both  the 
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phase  angle  W  of  the  eomplex  stress  intensity  factor  (K)'  and  the  ratio  of  stress  components  o.,  and  o„  along 
the  bond  line  are  invariant  with  respect  to  distance  from  the  crack  tip  an  are  re  a  e  . 
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Here  K  is  the  complex  stress  intensity  factor,  and  can  be  expressed  in  either 
K  =  Kj  +  K„  =  Ke''*').  This  expression  holds  in  the  near  tip  or  asymptotic  region 
singularity  eclipses  the  other  terms  in  the  Williams  expansion. 


rectangular  or  polar  form  (i.e., 
,  where  the  term  with  the  x'''' 


Also,  J,  the  contour  integral,  and  K,  the  magnitude  of  the  complex  stress  intensity  factor,  are  related  through 
an  effective  plane  strain  modulus: 
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Because  v  =  v,  =  1/2  for  incompressible  materials,  additional  simplifications  could  be  made.  More  general 
expressions  for  stress  and  displacement  components  in  the  near  tip  region  have  been  derived  ,  but  for  this  work 
only  eqs.  (1)  and  (2)  are  needed. 

>JTTMF.PTCAT.  MODELS 

A  typical  finite  element  model  is  shown  in  Fig.  2.  The  geometry  of  the  numerical  models  simulates  the 
experiiTntal  specimens  from  the  related  photoelastic  study.  The  specimens  are  glued  to  alummm  grips  which 
are  free  to  rotate  when  loaded.  Mode  mixity  is  varied  by  changing  the  crack  orientation  respert  to^the 
mode  I  loading  direction.  When  the  crack  is  oriented  horizontally,  mode  I  loading  tekes  place  and  K,,  . 

However,  as  the  crack  orientation  angle  a  varies,  the  mode  II  component  becomes  mcreasmgly  ' 

Crack  orientation  angles  of  0°,  15°,  30°,  and  45°  were  studied  in  both  the  experimental  and  numerical  work 

mode  minifies  could  be  studied  (die  corresponding  m«^^  about 

0°  to  30°).  The  numerical  models  also  allow  for  a  range  of  modulus  ratios:  E^/E,  -  1  (a  homogeneous 
specimen),  as  well  as  E2/Ej=  2,  4,  and  6. 

One  other  noteworthy  aspect  of  the  models  is  the  use  of  hybrid  elements.  The  use  of  incompressible 
materials  in  plane  strain  causes  an  indeterminacy  due  to  a  vanishing  of  a  term  (l-2v)  that  appears  in  the 
constitutive  equations.  This  results  in  a  singular  global  stiffness  matrix.  The  problem  is  resolved 
hybrid  formulation.  In  this  formulation,  the  hydrostatic  portion  of  the  stress  tensor  is  used  as  a  so  u  lo 
in  addition  to  the  components  of  nodal  displacements. 


•Throughout  the  text,  boldface  type  is  used  to  indicate  that  a  symbol  represents  either  a  complex  or  vector 
quantity. 


3 


j.pTFBMmmn  FRAr.TURE  papamf.tf.rs  from  numerical  daia 
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Mag  a  polynomial  carve  to  the  data  so  that  tan-(o./>)  «  a  function  of  r  Hie  constant  tem.  of  this 
polynomial  corresponds  to  t3n'\ajc^)  at  r  -  0,  and  is  an  estimate  o  . 

This  regression  method  has  been  used  previously  with  similar  numerical  models  and  comp^ed  wth 

::l«:  nL"™des  on  the  bond  line,  die  mean  of  the  tesute  for  each  element  is  determined  for 
use  with  the  regression  analysis. 

Equation  (2)  is  used  to  detennine  K,  the  magnimde  of  the  complex  stress  h,^ 

!  The  narameter  J  is  found  because  of  its  ease  of  evaluation  m  numerical  models  and  because  ot  its  roDust 

r  iH  vflnahle<;  of  the  elements  that  the  area  contains.  The  value  of  K  is  then  found  usmg  q.  C  )• 

wL  tli  K  has  also  indicated  good  agreement  with  existmg  experimenuii 

data 
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RESULTS 

Crack  orientations  of  0°,  15°,  30°,  and  45°  were  considered.  Etotic  \  ^2698  psi 

were  used  to  test  the  effects  of  relative  stiffness.  The  top  portions  of  the  specimens  always  have  E  f)98  ps 
(18  60  MPa),  corresponding  to  the  stiffness  of  a  photoelastic  polymer  used  m  related  ^23  4 

[Moduli  of  tl!;  lowerUon  of  the  specimens  were  varied  to  ch-ge  ^  ^ 

M)  were  applied  to  the  aluminum  grips,  givmg  a  nominal  stress  of  7.00  psi  (48.3  kP  ). 

Figure  3  shows  the  variation  in  K  with  E,/E,  for  all  crack  orientations. 
elevations  in  K  of  15-17%  occur.  However,  changes  in  Ej/E,  above  this  level  have  little  effect  o  .  y, 

Fig.  4  shows  the  variation  in  T,  the  phase  angle  of  K.  As  E,/E,  changes  from  to  2  J ^  J  ^ 
worst  case  (i.e.,  for  a  crack  orientation  of  0°).  Variations  in  the  phase  angle  are  smaller  as  E,/E,  changes  fro 

2  to  6. 

interface  due  to  diffusion  of  the  two  materials  during  the  rnanufacmring  process. 

The  primaiy  reason  for  the  weakness  of  the  E,/E,  effect  is  that  Ihe  modulus  mismatch  »  toge  changes 
o..  only,  with  only  small  changes  in  ^  »d  o  f  f „f  a  "  *  rimatch 

“nn  whS"" 

ill  the  more  compliant  material  and  compressive  in  the  *  ^„'from  a  cracked 

plots  o„  vs.  y  as  the  bond  line  is  traversed  for  a  specimen  with  EJE,  -  6.  The  data  is  taken  irom 


V 


,  ♦  frr.m  thp  ‘inecimen  boundaries  and  from  the  crack  tip. 

a  values  thus  change  substantially  with  Ej/E„  whereas  „ 

not  affect 'E  Changes  effected  in 'E  by  altering 

Since  'E  relates  to  0,/a,  in  the  near  the  near  tip  region.  A  similar  effect  occurs 

relative  stiffness  are  induced  by  sma  1  f  related  through  eq.  (2)).  Here,  a,,  is  a  fraction 

with  the  magnitude  K  or  with  the  valu  contribution  due  to  differential  contraction  being  a 

c„mpo„.m  =lo„g  the  concur  htet  '’J  'o^pS,  "Lmbuttons  to  different  pnm  of  the  contour, 

function  of  y.  This  fraction  component  makes  equal  o  w 
so  that  no  net  increase  in  J  is  produced. 


SUMMARY  AND  CONCLUSIONS 


cracks  along  interfaces  between  — i^ 

motors  and  have  a  bimaterial  parameter  ^  ^om  numerical  results  using  bond  line  fraction 

angle  of  the  complex  stress  intensity  factor  determined  using  J  integral  determinations 

data.  The  magnitude  of  the  complex  s^ess  mten  ^  magnitude  nor  the  phase  angle  varied  sipificantly 

For  a  variety  of  mode  mixities,  it  was  foun  ^  ^  substantial  changes  took  place  as 

as  the  modulus  ratio  for  the  material  coinbina  mbinations  of  propellant  and  liner  used  in  solid  rocket 

the  modulus  ratio  w^  je  of  kL  functions  of  loading,  geometry,  and  crack  length,  and  are  not 

motors  the  magnitude  and  phase  angle  o  n  ^  h  i  in  pr 

strongly  influenced  by  the  relative  stiffness  of  the  prope  an 


Young's  modulus 
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K 
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J  integral  value 

magnitude  of  the  complex  stress  intensity  factor 
complex  stress  intensity  factor 
bimaterial  parameter 


y  Poisson's  ratio 

fr,, .  <^xy .  ^yy  components 

xp  ■*  phase  angle  of  the  complex  stress  intensity  factor 
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Fig.  2  -  Typical  Finite  Element  Model  for  a  Bimaterial  Specimen 
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